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Abstract 

Previous work on dynamical black hole instability is further elucidated within the 
Hamilton-Jacobi method for horizon tunneling and the reconstruction of the classical 
action by means of the null-expansion method. Everything is based on two natural 
requirements, namely that the tunneling rate is an observable and therefore it must 
be based on invariantly defined quantities, and that coordinate systems which do not 
cover the horizon should not be admitted. These simple observations can help to 
clarify some ambiguities, like the doubling of the temperature occurring in the static 
case when using singular coordinates and the role, if any, of the temporal contribution 
of the action to the emission rate. The formalism is also applied to FRW cosmological 
models, where it is observed that it predicts the positivity of the temperature naturally, 
without further assumptions on the sign of energy. 
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1 Introduction 



In previous papers [H E] we considered the quantum instability of dynamical black hole using 
a variant of the tunneling method introduced by Parikh and Wilczek (PW) in the static case 
[3|,|1] to uncover aspects of the back reaction effects; see also [I5J. The PW tunneling method 
was refined and extended to more general cases in [6] and other papers [H [8] , including its 
relation with thermodynamics [9]. Other points of view and criticisms can be found in [10J. 

In the present paper, namely in section [2} we review this variant method, which in the 
static case goes under the name of Hamilton- Jacobi (HJ) method of tunneling, and was 
introduced in [TTJ [12]. A complete comparison analysis with the Parikh- Wilczek method 
has been presented in [13J. The HJ tunneling (across a horizon) connects the invariant 
surface gravity to a dynamical local temperature through the leading term in the black hole 
tunneling rate ( 11.11) . 

First, we would like to notice that the main novelty of the HJ method is its manifest 
covariance, compared to the original Parikh- Wilzcek approach [31 H] where a non manifestly 
covariant Hamiltonian formulation was used. With regard to this issue, it is our opinion 
that HJ method is particularly suitable for the generalization to the dynamical spherically 
symmetric case, as we will try to substantiate in this paper. 

The tunneling method provides not only new physical insight to an understanding of the 
black hole radiation, but is also a powerful and simple way to arrive at an expression for 
the surface gravity for a vast range of solutions, especially non-stationary and, in our case, 
dynamical black holes. This looks important since several definitions of the surface gravity 
for evolving horizons have been proposed in the past, all fitting some kind of first law of black 
hole mechanics more or less equally well. Still different results are advocated for expanding 
cosmological black holes in [TH [15] . A comparison is discussed thoroughly in [161 E] > but 
we anticipate that Hayward's definition [T5] of surface gravity, related to Kodama's theory 
[T9] of spherically symmetric space-times, is the most interesting to us since it is the one to 
which the tunneling method naturally leads. 

As we said above, a way to understand Hawking radiation is by means of tunneling of 
particles through black-hole horizons. Such a tunneling approach uses the fact that the 
WKB approximation of the tunneling probability for the classically forbidden trajectory 
from inside to outside the horizon is: 



where / is the classical action of the (massless) particle, to leading order in h. It is of the 
utmost importance that the exponent be a scalar invariant, otherwise no physical mean- 
ing can be given to T. If, in particular, it has the form of a thermal emission spectrum 
with 21m / = /3u, then the inverse temperature (3 and the particle's energy uj have to be 
separately scalars, since otherwise no invariant meaning can be given to the horizon tem- 
perature, which would not be an observable. If, even in the presence of reasonable physical 
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conditions, more than one prescription for defining an invariant energy is available, then also 
more notions of invariant temperature will exist and further analysis, or observations, could 
be needed. In most cases this could only be a scale transformation or a choice of a different 
family of observers. 

In principle, all the Standard Model particles are expected in the Hawking radiation spec- 
trum. However, most of the calculations in the literature have been performed just for scalar 
fields. Spin one-half emission was considered in [20] for stationary black holes, in [21] for the 
special case of the BTZ black hole and finally [22] studied the case of evolving horizons. 
In Ref. [2] we have already discussed the Hamilton- Jacobi tunneling method for arbitrary 
spherically symmetric dynamical black holes. In this paper, we would like to present a sys- 
tematic derivation of the results, pointing out the key points and making use of a covariant 
coordinate approach^ which is particularly convenient in the discussion of several explicit 
examples we will consider as applications of the general method. We point out that our 
description of the tunneling rate is not equivalent to the "canonically invariant" description 
due to Chowdhury [24], which in fact obtains half the correct temperature for black holes 
and infalling shells. Probably a comparison analysis would be worthy, at this point. 

Recall that any spherically symmetric metric can locally be expressed in the form 

ds 2 = % j (x i )dx l dx i + R 2 (x i )dQ 2 , i,j e {0, 1} , (1.2) 

where the two-dimensional metric 

d^ = ltj (x i )dx i dx j (1.3) 

is referred to as the normal metric, x l are associated coordinates and R(x l ) is the areal 
radius, considered as a scalar field in the normal two-dimensional space. Another relevant 
scalar quantity on this normal space is 

x(x) = i ij (x)d i Rd j R. (1.4) 

The dynamical trapping horizon, say H, may be defined by 

X (x) =0, d lX ^0. (1.5) 



H 



H 



This is equivalent to the vanishing of the expansion 9 m of the null, future directed congruence 
which is normal to a section of the horizon. In general there will be another null congruence 
of "incoming rays" with a related expansion 9/ n y. if 0( n ) < and also C n 9^ < along H, 
then the horizon is of the future, outer type. Most of our results will be valid for this type 
of trapping horizon. 

The Misner-Sharp gravitational mass, in units G = 1, is defined by 

m(x) = ^R(x) (1 - x(x)) ■ (1.6) 



1 For static black holes a coordinate- free formulation has been given by Stotyn and co-workers 
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This is an invariant quantity on the normal space. Note also that, on the horizon, m|# = 
m = Rh/2. Furthermore, one can introduce a dynamic surface gravity [2] associated with 
this dynamical horizon, given by the normal-space scalar 



H 



di^^djR) . (1.7) 



H 



Recall that, in the spherical symmetric dynamical case, it is possible to introduce the Kodama 
vector field K, such that (K a G a/3 y (S = that can be taken as its defining property. Given 
the metric (11.21) . the Kodama vector components are 

K l (x) = -^=^djR , K e = = . (1.8) 

The Kodama vector gives a preferred flow of time and in this sense it generalizes the flow 
of time given by the static Killing vector in the static case. As a consequence, we may 
introduce the invariant energy associated with a particle by means of the scalar quantity on 
the normal space 

u = -K i d i I ) (1.9) 

where I is the classical action of the massless particle, which we assume to satisfy the 
Hamilton- Jacobi equation 

Y j dildjl = 0. (1. 10) 

We stressed above the importance to have at disposal an invariant definition of energy. 
Eq. (11.91) certainly satisfies this requirement if the action is, as it normally is, a scalar. In the 
following our aim will be not so much a detailed picture of the physical process of horizon 
tunneling, as to show that there is a precise invariant prescription to deal with the imagi- 
nary part of the action, in case there is one, which is valid for all solutions in all coordinates 
systems which are regular across the horizon. 

The plan of the paper is the following. In Section II we discuss the role of coordinate invari- 
ance in the Hamilton- Jacobi method. In Section III, the first law for dynamical black holes 
is derived. In Section IV, FRW space-times are investigated in detail as relevant examples 
of dynamical horizons. The naturally obtained positivity of the temperature is emphasized. 
In Section V, static black holes are revisited, and the issue of Hayward surface gravity and 
Killing surface gravity is discussed in detail. The paper ends with the Conclusions. 



2 Hamilton-Jacobi Tunneling and Gauge Invariance 

To begin with, we would like to point out that, within this context, the tunneling has to 
be interpreted in a different way from the usual textbook treatment of tunneling effect in 
(non-relativistic) quantum mechanics. In fact, strictly speaking, there is no barrier here: 
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an issue which has been clarified by Parikh and dubbed by him as "secret tunneling" [1]. 
Furthermore, the appearance of an imaginary part in the "classical" action is due to the 
presence of a dynamical horizon and the use of Feynman's ie-prescription in dealing with 
an otherwise divergent integral within the classical realm. This is the way, together with 
equation ( 11.1)1 . in which quantum effects enter into the game. 

To be more specific, our key assumption is the following: that we can reconstruct the 
whole action J, from <9;7 by means of 

dx^J, (II. 1) 

where 7 is an oriented, null curve, with at least one point on the dynamical horizon. We 
can split the integration along this null curve in two pieces, one very near to the dynamical 
horizon, the remaining contribution living in the regular part of space-time. Then, we may 
perform a near-horizon approximation in the first integral. As a consequence, we have to use 
a system of coordinates which is regular on the horizon, otherwise this procedure cannot be 
applied, due to highly singular quantities involved. This procedure permits to evaluate the 
integration along what we can call, certainly with some abuse of notation, the "radial" and 
"temporal" parts and see if an imaginary part shows up, and what is their relation. 

In Ref. [2], we have presented a derivation of the relevant imaginary part of the classical 
action and tunneling rate ( II. lj) . valid for a future trapped horizon, which reads 

Im / = Im {^J ' dx%rj = ^ , (II.2) 

where kh is the dynamical surface gravity ( 11.71) . and ujh the Kodama energy evaluated on 
the dynamical horizon. These quantities are scalars in the normal space, thus the leading 
term of the tunneling rate is invariant, as we expect to be from an observable. 
For the sake of completeness and clarity, we are briefly reporting the essential steps of the 
derivation in several coordinate systems. 



The EFB Gauge. Let us start with giving the computation as in [2]. The key point is the 
observation that it is always possible to rewrite, locally, any spherically symmetric metric in 
an Eddington-Finkelstein-Bardeen (EFB) form which is regular on a trapping horizon. For 
black holes, one is concerned with future trapping horizons and so the advanced (rather than 
retarded) form is used: 

ds 2 = -e^Cdv 2 + 2e*dvdr + r W , (II.3) 

where we are using x l = (v, r) as coordinates, and C — C(v, r), ^ = ^(v , r). In this gauge, 
life is easy since the areal radius of the spheres of symmetry R turns out to be equal to 
the r-coordinate, i.e. R = r, and x = C: the horizon location is defined by the condition 
C(v,t)\h = 0. The Kodama vector and the invariant energy assume the simple expressions 
given by K = (e - *, ), and uj = —e~^d v I while the invariant surface gravity is just given 
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by Kh — dj.Cn/ 2- One may note that \1/ transforms as an ordinary Liouville field, i.e. 
^ — > ^ + In \ dv/dv |, under v — > f5(v), making invariant under reparametrizations of the 
advanced time coordinate. 

Expanding along a null direction from a point on H outward in a neighbourhood of the 
horizon gives (since Ch = 0) 

= e* H ArAv, (II.4) 

which shows that the temporal (v-) contribution does not play any role in what concerns us 

at the moment (the evaluation of Im I) . 

From the HJ equation, we see that, for outgoing modes, 

C(d r I) = 2u . (II.5) 

Thus, taking (111.4ft and (III.5j) into account, one has 



ImJ = Im / (d r Idr + d v Idv) 

J -V 



[in / dr — (II.6) 
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2 Im / dr (II. 7) 

J 1 o r C\ H [r - r H - zO) 



(II.8) 



where again, the quantity C has been expanded around the horizon along the null direction, 
that is 

C{v,r)&d r C Ar + ... (II.9) 

and Feynman's ze-prescription has been implemented in order to deal with the simple pole 
(III.7p . kh = d r C\ H /2 coincides with our geometrical expectations and we see that, in EFB 
coordinate system, the temporal integration does not give any contribution to the imaginary 
part of the action of particles tunneling through the (trapping) horizon. 

However, for practical reasons, it might be convenient to work in other (regular on the 
horizon) coordinate systems. So, denoting the temporal and spatial coordinates by x l = 
(t,r), we are going to discuss three more instructive gauges. 

The r-Gauge. The normal metric here is non-diagonal, but as in EFB gauge, R = r. We 
have 

d s 2 = d~i 2 + r 2 dQ 2 , (11.10) 
where the reduced normal metric is 

dj 2 = -E(r,t)dt 2 + 2F(r,t)dtdr + G(r,t)dr 2 , F^O. (11.11) 
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The horizon is located where 

X (t,r)=Y 3 d l Rd J R = 1 rr (t,r) 



E 



EG + F 2 

vanishes, i.e. at Eh = 0, provided Fh ^ 0. The Kodama vector reads 
1 



K 



.VW+EG' 
and the invariant energy 
d t I 







U) 



VF 2 + EG 

The dynamical surface gravity is, from (11.71) . 



2F 3 



E'F - ^EG 



;iI.12) 



;iI.13) 



(11.14) 



;n.i5) 



1-1 



where an overdot and a prime denote differentiation with respect to t and r, respectively. 

Ar (the other solution being 



H 



From the metric the null, radial, expansion gives At = — ^ 

related to the ingoing null ray), so within the near- horizon approximation and after some 
calculation we get x — 2«ij( r — t h)] also, dtl = —Fhuj from definition (1II.14P and the horizon 
condition E = 0. 

Splitting the integration along 7 according to what we said above, we end with I given by the 
sum of a real term and a possibly imaginary part coming from near horizon approximation: 



/ I (drd r I + dtd t I) = I dr 

J 'v J <y 



d r I + -G h ujh 



(11.16) 



What is remarkable is that in this gauge, the temporal, t-part is present^], but being regular, 
it does not contribute to the imaginary part of the action. Making use of the HJ equation, 
the Kodama energy expression fill. 141) and equation fill. 121) as well, one has 



uF 



VEG + F 2 
thus, for outgoing modes, 
uF 



d r I - uj 2 G = 



d r I 



VEG + F 2 X 



(2 + 0(x)) • 



;n.i7) 



(11.18) 



2 Since it contributes to the total action through the \{Goj)h term. 
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Making use of this equation and Feynman's prescription, Eh = and \ — 2^n{ r — r n), one 
has for the outgoing mode 



Im / = Im / dr d r I 
'i 



= Im £ 



uF 1 + y/l + 0( X ) 1 

1m / dr 



'7 
7TU H 



VF 2 + EG 2 Kfj (r-r H - iO) 

(11.19) 



in agreement with the EFB gauge. 

The Synchronous Gauge. The second coordinate system we would like to consider is 
described by the line element 

ds 2 = -dt 2 + -J-^dr 2 + R 2 {r, t)dtt 2 = d^ 2 + R 2 {r, t)dQ 2 , (11.20) 
±> [r, t) 

in which the metric is diagonal, but the proper radius of the spheres of symmetry R is a 
function of the coordinates r and t. In this case, one has 

X = -(dtR) 2 + B(d r R) 2 , (11.21) 

thus the future sheet of the trapping horizon xh = is given by 

(d t R) H = -y^(d r R) H , (11.22) 

in which we are assuming again a regular coordinate system on the horizon, namely that Bh 
and its partial derivatives are non-vanishing. The Kodama vector reads 



K = (V Bd r R, —VBd t R, 0, 0) , (11.23) 
and the invariant energy 

u = y/B(-d r RdtI + d t Rd r I) . (11.24) 
The dynamical surface gravity may be evaluated and reads 

kh = \ (-1d 2 t R H + 2B H d 2 R H + ±-d t R H d t B H + d r R H d r B H ^j . (11.25) 

Making use of the horizon condition, we may rewrite 

k h = \ (-2d 2 R H + 2B H d 2 R H - -±=d r R H d t B H + d r R H d r B H ^j . (11.26) 

From the metric, the HJ equation reads simply 

- {d t lf + B{d r I) 2 = 0. (11.27) 
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As a consequence, the outgoing temporal contribution is equal to the radial one and we have 
1 = 2 f drd r I . (11.28) 

The HJ equation and the expression for the invariant energy lead to 

8 r I = . (11.29) 

Bd r R + yfBd t R 

which has a pole on H . Making the expansion along the outgoing null curve, for which 
At = — Jg- Ar, in the near-horizon approximation, one gets 

Im I = 2 • Im I dr - = — , (11.30) 

] 1 2n H {r - r H - i0) k h 

which again coincides with the previous result. But notice that in this gauge, the temporal 
contribution is essential indeed to provide the correct result: without it the temperature 
would be doubled. 

Conformal 2D Gauge. Another coordinate system where the temporal contribution to 
the action plays an essential role is the general diagonal form of a spherically symmetric 
metric, which reads 

ds 2 = e^' r) {-dt 2 + dr 2 ) + R 2 (t, r)dn 2 . (11.31) 

In this form, the normal metric is conformally related to the 2-dimensional Minkowski space- 
time. The x function simply reads 

X = e~* ((-d t R) 2 + (d r R) 2 ) , (11.32) 
which leads to the (future trapped) horizon condition 

{d t R) H = -{d r R) H . (11.33) 
The Kodama vector and associated invariant energy are 

K = (d r R, -dtR, 0, 0) , (11.34) 

uj = e~* {-d r Rd t I + d t Rd r I) . (11.35) 
The dynamical surface gravity reads 

k h = \e-^ H {-d 2 t R + d 2 R) . (11.36) 

Due to conformal invariance, the HJ equation is the same as in two-dimensional Minkoswki 
space-time and we may take 

d+I = d t I + d r I = . (11.37) 
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Since the null expansion condition leads to Ax + = At + Ar = 0, for outgoing modes we get 



(drd T I + dtd t I) = 2 / drd r I . 



Furthermore, due to flll.35j) and flll.37j) . one has 



d r I 



e-^(d r R + d t R) 



(11.38) 



;il.39) 



and we have a pole at the horizon. Making use of near horizon approximation along the null 
direction, from (fIL33l) and fjlL36|) . one has (d r R) H + (d t R) H = 0, At + Ar = 0, thus 



d r R + d t R 



{dlR • 
2k h (t 



■r H ) + 



d 2 tt R + dlR) 



(11.40) 



As a result, making use of Feynman's prescription, one again arrives at equation ( 111.21) . 
In the previous computations various choices of signs have been applied in such a way that 
it may seem they were chosen somewhat ad hoc in order to get the wanted result. This 
is not so. Once the future sheet of the trapping horizon has been chosen, and the sign of 
the Kodama vector so determined that it is future directed, no other sign uncertainties will 
occur for either outgoing or ingoing particles. On the other hand, if there exist a past sheet 
in the trapping horizon then using the tunneling picture we may as well compute the action 
along an inward directed^) curve at the horizon. Then there will be again a non-vanishing 
imaginary part, but we can interpret it as a small absorption probability. 



3 The First Law for Dynamical Black Holes 

Here we present, for the sake of completeness, a derivation of a version of the first law. To 
this aim, let us introduce another invariant in the normal space, related to the stress-energy 
tensor: 

T (2) = 7 *%. (III.l) 
First, let us prove the following invariant relation 

k h = 7^- + ZkRhTP , (IIL2) 
2K H 

valid on the dynamical horizon. 

We may use the EFB gauge, in which R = r and x = C and, in this gauge, we have 

K H = \drC H - (IIL3) 



3 The ambiguity inherent in this and analogous terms is easily resolved if the manifold is asymptotically 
flat. 
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On the other hand, the Einstein equations are very simple in this gauge and we have 

' ~ ° - -0 r C = -A7ir 2 T v v . (III.4) 



2 2 
Thus, on the horizon we get 

\ ~ r ~fd r C H = -^r\Tf (III.5) 
since T r r H = 0, which leads immediately to flIII.2p . 

However, it is easy to prove the same relation in another coordinate system, for example the 
r-gauge. To this purpose, let us introduce the horizon area and the areal volume associated 
with the horizon, with their respective differentials: 

Ah = 47rR H , dA H = 8nR H dR H , (III.6) 

V H = -7iR H , dV H = 4nR 2 H dR H . (III. 7) 

3 

Then a direct calculation gives 



^dA H = d \-^\ + T%>dV H . (III.8) 

In turn, this equation can be recast in the form of a first law, once we introduce the MS 
energy at the horizon: 



dm = —d {-f - T^dVu , (III.9) 
where Sh = Ah/ '4 generalizes the Bekenstein-Hawking black hole entropy . 

4 An Explicit Example: The FRW Spacetime 

A very interesting example of the tunneling method is provided by a generic FRW space-time 
with flat spatial sections, 

ds 2 = -dt 2 + a(t) 2 dr 2 + a(t) 2 r 2 dtt 2 . (IV.l) 

At first glance, this example can seem to lie somehow outside the main stream of the paper 
which, up to this point, has been devoted to the study of dynamical black hole horizons. 
However, what we are considering is a truly dynamical horizon but, this time, of cosmological 
interest . 

As explained above, metric (lIV.ip belongs to the class of Lemaitre-Rylov, a subset of the 
synchronous dynamical space-times, and we have just showed this to be a coordinate system 
where the evaluation of the imaginary part of the action gets a contribution due to the 
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integration along the t-coordinate. 

The normal reduced metric is diagonal with coefficients 7^ = diag(— 1, a{t) 2 ) and % = 1 — 
r 2 d(t) 2 . The dynamical horizon is the Hubble horizon (we assume H(t) > 0, see also |25j ) 

r » = (IV - 2) 

with the Hubble parameter H(t) = a/ a. The dynamical surface gravity is 

k h = -\H + ^-\ , (IV.3) 




and the minus sign refers to the fact the horizon in question (1IV.2I) is, in Hayward's termi- 
nology, of the inner type. For example, in the Einstein-de Sitter model, with a{t) oc t 2 / 3 , we 
would obtain 



1 

k h = 



6f 



while the de Sitter model in inflationary coordinate has H = and kh = — a/A/3, though 
this result is valid in any other patch; furthermore, Kh — is only possible in a radiation 
dominated universe, where a(t) oc \/t. One easily sees that in general, for a flat model with 
a{t) oc t n one has 

so only for n < 1/2 is our surface gravity positive. This regime should not be physically 
allowed, however, since radiation dominated models occur either with massless particles or 
with ultra-relativistic massive ones, both of which are limiting cases. It seems as if in these 
cases we should define the temperature as T = \kh\/2tv. However, we will see that the 
tunneling method just gives the right signs without invoking absolute values. 

For a massless particle, the reduced HJ equation is 

d t I = ^-d r I. (IV .4) 

a{t) 

The full classical action of outgoing particles is 

dx i 0J , (IV.5) 
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with 7 an oriented curve with positive orientation along the increasing values of x % = (t, r). 
Radially moving massless particles follow, of course, a null direction. Then, we can perform 
a null-horizon radial expansion, 

= ds 2 = -dt 2 + a 2 (t)dr 2 (IV.6) 
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which gives 

At = a(t)Ar. 



(IV.7) 



for out particles. The outgoing particle action, that is the action for particles coming out of 
the horizon, is then 



dtd t I+ / drd r I 



2 / drdj 



(IV.8) 
(IV.9) 



upon using equation (IIV.7I) . The Kodama vector reads K — (1, — r77, 0, 0) and the invariant 
energy of a particle is given by 



u = -3J + rHd r I . 



Thus 



dj 



a(t)uj 



ra(t)H-l 
and we have that 

J, = -2 / dr 



a{t)uj 



1 - ra{t)H{t) 



(IV. 10) 



(IV.11) 



(IV. 12) 



Expanding the function f(r,t) := 1 — ra(t)H(t) close to the horizon, again along a null 
direction, we get 



/M 



77 



2a(t) 



H 



277 



[r - r H ) + 



= 2K H a(t)(r - r H ) 



(IV. 13) 



where kh represents the (dynamical) surface gravity associated to the horizon. 
The action of the outgoing particle now reads, 



dr 



K H {r -r H - i0)' 



(IV. 14) 



In order to deal with the simple pole in the integrand, we implement the Feynman's ie - pre- 
scription. In the final result, beside a real (irrelevant) contribution, we obtain the following 
imaginary part: 



ImL 



(IV.15) 
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As a consequence, we may interpret T = — Kh/2tt > as the dynamical temperature 
associated with FRW flat space-times. In particular, this gives naturally a positive tempera- 
ture for de Sitter spacetime, a long debated question years ago, usually resolved by changing 
the sign of the horizon's energy. It should be noted that in the literature, the dynamical 
temperature is usually given in the form T = H/2tt, with a missing term depending on H, 
exceptions being the papers [26]. Again, T becomes negative only for the unphysical flat 
models with n < 1/2, or perhaps we may say there can be no tunneling processes from them. 



It is instructive to reconsider the FRW tunneling computation in another coordinate system 
discussed in previous Sections. Making the coordinate change R := ra(t), the metric assumes 
the form of r- gauge, namely 

ds 2 = -(1 - H 2 R 2 )dt 2 - 2HRdtdR + dR 2 + R 2 dtt 2 (IV.16) 

Note that the metric remains regular on the horizon and the associated normal metric is of 
the type fill. 111) . As a result, \ = j RR and the dynamical horizon is defined by H 2 R 2 H = 1. 
Of course, the dynamical surface gravity has to remain unchanged with respect to (liV3j) . 



However, in this particular case, the Kodama vector is very simple, K = (1, ). As a 
consequence, the invariant energy is just uo = —d t I, and the HJ equation for a massless 
particle along a radial trajectory reads: 

- (d t I) 2 - 2HR(d t I)(d R I) + (1 - H 2 R 2 )(d R I) 2 = . (IV. 17) 

We know already from §11, that the integration along temporal coordinates gives merely a 
real contribution. Thus, an imaginary contribution to the particle action, comes only from 
integration along the radial direction. The HJ equation (II V. 171) supplemented of the Kodama 
energy constraint, gives 



(d R ih - r^jpR- (IV - 18) 

Making a null-expansion on the horizon, 

= As 2 = -2AtAR + AR 2 (IV. 19) 
we see that 

(1 - H(t) 2 R 2 ) ^-2\H+^—\AR = 2k h (R - R H ). (IV.20) 




We finally get the expected results, 

f u,Hr(i + y/1 + 0(R-R H j) 

lmL=- dR ^ L = ^^2L (iv.21) 

J -2(H + H/2H)(R-R H -iO) k h 

with kh provided by (IIV.20I) . 
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Everything we said above generalizes straightforwardly to models with non- vanishing spatial 
curvature (k = 0, ±1), except that the surface gravity is given by the more complicated 
formula 

kh = - (h 2 + ±H+J^j R H (IV.22) 

where Rh = (H 2 + k/a 2 )~ 1 ^ 2 is the cosmo logical trapping horizon, coinciding with the Hubble 
sphere in the flat case. In this case, one may note that k h = is possible as soon as 

a(t) = \l-kt 2 + cit + c 2 , (IV.23) 

where c\ and c 2 are constants. Of course, for k = +1, the solution is real only in a finite 
range between a big bang and a big crunch. 



5 Static Black Holes 

Static black hole solutions may be considered as a special case of dynamical ones. However, 
they are consistent solutions of relativistic theory, Einstein or modified alternative theories, 
and they deserve a separate analysis, even though the horizon tunneling, which we are going 
to discuss according to the general procedure outlined in the previous Sections, is a signal 
of their quantum instability and the static hypothesis we assume is only an approximation, 
strictly valid only for limited periods of time. 



5.1 Kodama-Hayward vs Killing Surface Gravity 

Let us consider the Schwarzschild diagonal gauge, which can be written as 

dv^ 

ds 2 = -V(r)dt 2 + — — + rW . (V.l) 
W(r) 

In this case the function x defined in (11.40 coincides with W. Thus, if we assume that V and 
W have the same simple zeros, we have, on the horizon defined by \ = W = 0, 

Vff V'rr 

w H -w H < v - 2 > 

via de l'Hospital ratio rule. These coordinates, as is well known, are singular on the horizon, 
so that the null expansion we have to use is meaningless, the temporal contribution to the 
action being ill defined. The use of these singular coordinates has been the origin of a large 
number of papers, containing several proposals to deal with the ambiguity [10] • Recall that 
in the original paper by Parikh & Wilczek on the tunneling method [3], a clever use of a 
coordinate system, regular on the horizon, known as Painleve-Gullstrand coordinates (PG), 
was advocated. The general PG gauge reads 



r/.s 2 = - \ 'J1- - 2 \/ ^(1 - W)drdt + dr 2 + r 2 dVi 2 . (V.3) 
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According to our assumption (1V.2I) . this gauge is indeed regular on the horizon, being of 
what we termed r-gauge type. A variant of this gauge is the EF gauge, where an advanced 
(retarded) time appears, i.e. 

ds 2 = -Vdt 2 -2\ —drdt + r 2 dn 2 . (V.4) 

v w 

In both gauges, we may apply the general analysis and conclude that the temporal contri- 
te 



bution is absent. The Kodama vector is K = (y-nr, ) and the related invariant energy 



to = y ^yd t I. One should note that the Kodama vector and the invariant energy do not 
coincide with the Killing vector and dtl unless V = W. This has some consequences, since 
the general theory gives a surface gravity of 

«* = , (V.5) 



instead of the surface gravity associated to the Killing vector <9 ( 



(V.6) 

More details can be found in [22], but for now we just limit ourselves to note that the 
tunneling probability 

r ~ e -2 "^ , (V.7) 
is the measurable quantity and since 

^ = -, (V.8) 

where E = dtl is the Killing energy, no contradiction within the static approximation can 
be found. However, as soon as the space-time becomes dynamic, a Killing vector is useless 
and the only invariant energy is u, indeed. 

5.2 Dilaton-Maxwell-Einstein Black holes 

In this Subsection, we would like to give a brief review to black hole static solutions in the 
form (jVTjl . with W ^ V. 



As far as one is concerned with Einstein's GR, solutions of this kind are forbidden by well- 
known uniqueness theorems, cfr. [2B] ■ However, it is not difficult to face them in alternative 
theories of gravity, a label which actually contains a lot of different material. In order to be 
more precise, let us consider the specific example provided by the so-called Dilaton-Maxwell- 
Gravity (DMG) as in [13] and [29, [301 EH E2]- We typically start by an action such as 



J d 4 x^ (R - 2(V0) 2 - V(<f)) - e~ 2 ^F 2 ) 



(V.9) 
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where R is the scalar curvature, F 2 = F a i J F ab and £ governs the coupling of the dilaton 
with the Maxwell field. Varying the action (IV. 91) with respect to the metric, Maxwell and 
dilaton fields yields the EOM for the respective fields. It is easier if we consider Maxwell 
fields generated by an isolated electric charge. Then, the ansatz 

ds 2 = -U{r)dt 2 + U-\r)dr 2 + H 2 {r)dtt 2 (V.10) 

for the line element, satisfies the EOM for the metric under very general conditions]^ 
Assuming b < a, a first class of solutions is given by [30J 

1 - - 

U(r) = 1 & H(r)=r 2 -br, (V.ll) 

r 

another one is instead 

U(r) = l--, & H(r) = r 2 -br. (V.12) 
r 

For a quite general class of solutions, but without any attempt of universality, the U(r) 
function in (IV. 101) can be expressed in terms of H(r) and the dilaton as 

U ^ = H 2 (r)4>'(r)+^H(r)H'(r)' 

where rj is an integration constant and the prime denotes derivation with respect to the 
argument. 

Next, consider a conformal transformation 

ds 2 = n 2 ds 2 (V.14) 
so that 

= J^ry ( v - 15 ) 

The conformal metric looks like (IV. II) . 

ds 2 = -V(r)dt 2 + W^irjdr 2 + r 2 duj\ (V.16) 

with 

= r 2 U(r) = r 2 ^r + V ) 
y ' H 2 {r) H*tu + £H S H> y J 

W(r) = mim = J^±±_ (V18) 
If the spacetime possesses a horizon, this will be located at r , s.t. W(r ) = 0, i.e. where 



4 The Maxwell field then reads F tr = ^ , Q being the electric charge. 
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1. U(r ) = 0, that is r = 



-v/t > 0; 



2. H{tq) = 0. Note that in this case 

lim V(r) = oo lim Q 2 (r) = oo , (V.19) 

r— >ro r— »ro 



that is, the conformal transformation becomes singular. 

It is simply a question of algebra to compute the Kodama-Hayward surface gravity and see 
that it is well-defined in both cases 1 and 2. But, with regard to the Killing surface gravity, 
we have to distinguish carefully between the two cases under examination. Indeed, case 1 is 
easy to treat and gives k,k = U' /2: something we had to expect from the very beginning in 
consideration of the well-established result according to which the Killing Hawking temper- 
ature = kk/2tt is conformally invariant |33j. But things go radically different in case 2 of 
a singular conformal transformation, where the same definition of a Killing surface gravity, 
positive and finite on the horizon becomes, at least in the general case, questionable. 

This analysis, however, is sufficient to shed new light on the stringy black hole puzzle. 
Start in fact from the metric (IV. 101) with, for example, U(r) — 1 — a/r, a — const and 
H(r) = v/r(r — b), b < a. Perform a conformal transformation in order to get the GHS 
solution (30]. Since b < a, r = a is still an event horizon. The Killing surface gravity, 
being invariant under conformal transformation, does not feel the new physics introduced 
by the conformal transformation. Thus no change has to be expected in the extremal limit. 
The story is different for Hayward's surface gravity that goes like kh oc Hq and vanishes 
whenever the conformal factor vanishes, e.g. in the extremal limit for GHS solution. 



5.3 The Lemaitre-Rylov Gauge 

As further example, let us consider the Schwarzschild spacetime in coordinates (t, r, 6, ip) 
such that the line element can be expressed as 

ds 2 = -dt 2 + + (r g B) 2 dQ 2 , (V.20) 
where r g = 2m is the usual gravitational radius, and 



(V.21) 

We shall refer to these coordinates as the Lemaitre-Rylov gauge. This is indeed an interesting 
(time-dependent) gauge since - contrary for example to isotropic coordinates - (t, r) extend 
beyond the gravitational radius, r < r g . It is an explicit example of synchronous gauge. 
Notice further that [28J: 

1. The spacetime singularity is located at r = t; 



B(t,r) 



— (r - 1) 
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2. The horizon is located in correspondence of 

2r 

B(t,r)\ H = l, (r H -t) = -f. (V.22) 

3. Outgoing particles are such that % < 0; ingoing particles instead have % > (cfr. 
spacetime diagram in [28J for example). 

A detailed calculation can now be provided. 
According to the general procedure outlined above, 

7«=f" 1 b-i) , v^=5- 1/2 , (V.23) 



and 



r = ( ^ ° ) • (V.24) 

Making use of (II. 7j) . a direct computation leads to /c# = as expected. The explicit form 
of the Kodama vector K in terms of the Lemaitre gauge is 

#=(1,1,0,0), (V.25) 
so that the particle's energy is 

u = -d t I - d r I . (V.26) 
The H-J equation for radially moving particles is 

-(d t I) 2 + B(d r I) 2 = 0, (V.27) 
and for the outgoing particles, we take 

d r I = — ^£ . (V.28) 
y/B 



Thus 

9rI = ^±H. (V.29) 

£5 — 1 

The null expantion condition gives for the outgoing particle At = — namely we have 



/ = 2 / dr(O r I) + = 2 I (//w V ^_ + 1 1 - (V.30) 
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Let us define £/(t,r) := B — 1 and expand it along the null, radial, outgoing geodesic close 
to the horizon: 

B - 1 ~ d r B H Ar + d t B H At = 2d r B H Ar = — (r - r H ) . (V.31) 



r 9 



As a result 



I = r„ ( — 1) =^ ImI = ™ = AnmuJ /y_32) 
Jy{r-r H - 10) fe H 

which provides the correct Hawking temperature from the hole. The Lemaitre-Rylov gauge 
can be generalized beyond Schwarzschild spacetime. In general the metric we shall deal with 
is 

rfs 2 = -eft 2 + (l - F(i2))dr 2 + flW , (V.33) 
where V(R) and the areal radius R are function of r ± £, via the inversion of the relation 

r±t = G(R) = TjdR^(j^-VT=Vy (V.34) 

We are dealing with a dynamical synchronous gauge and we may apply the general formalism 
seen above. A straightforward calculation leads again to (jVJjpE. 



6 Conclusions 

In the last few years, many different proposals have been suggested in order to give a uni- 
versal prescription for the Hawking temperature of certain dynamical spacetimes, especially 
those endowed with future trapping horizons. In such cases, for example, one lacks the 
Kubo-Martin-Schwinger (KMS) condition so successful for equilibrium states. Nor is there 
generally available a analytic continuation to Euclidean signature within which one can judge 
the periodicity of Euclidean time, essentially equivalent to the KMS condition. In fact, one 
can even doubt that a temperature with the usual meaning is generally possible or useful 
at all. The final result is that time by time many prescriptions were proposed which either 
were not applicable to all the desired cases one can have in mind, or simply were invented 
ad hoc to keep the peace with the most well known cases. For example, for a slowly varying 
Schwarzschild black hole it seemed natural to keep the temperature equal to the instanta- 
neous value l/8irM(t) characterizing the static black hole, although this is not so obvious. 

On the other hand we are all accustomed to the remarkably universal properties exhibited 
by black holes so, moved by the wish of extending some of these to a dynamical regime, we 

5 As an exercise, the interested reader could be curious to see what happens by introducing Kruskal- 
Szekeres generalized coordinates for the metric (|V.1|) . The computation is rather long and tedius, but still 
confirms, as it must be, the general covariance of Hayward's surface gravity. 
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have written this note in order to clarify (hopefully) the so-called HJ method for horizon 
tunneling. 

The playground has been given by the class of spherically symmetric spacetimes, either 
static or dynamical, and within this class we can draw the conclusion that in order for 
the HJ method to work properly: (i) regular coordinates at the horizon are necessary; (ii) 
an appropriate notion of particle energy, that it should be a scalar, can be implemented 
at the level of HJ equation in such a way that imaginary contributions to particle action of 
outgoing particles arise after a near-horizon approximation; (iii) we have shown the relevance 
of the null-expansion of (to all purposes, massless) particles through the horizon in order 
to reconstruct fully the particle action. In turn, it is the particle action responsible for the 
particle production rate (11.11) . 

As a result, it should also be clear that a contribution to the rate from the temporal part 
of the integration, that is from the piece J d t Idt of the basic tunneling rate, is generally 
present but gauge dependent; it is only the full integration along "7* that produces a gauge 
invariant result. So for example, in certain gauges there is a temporal contribution, in some 
other gauge there is not. 

Finally, in static spacetimes violating the weak energy condition it seems that the choice 
between Killing and Kodama- Hay ward cannot be decided at the formal level, although it 
may be remarked that the difference amounts only to a trivial scaling, albeit one which 
can have less trivial effects on extremal black holes. It may also be noted that, as soon 
as a static black hole starts evaporating, the spacetime ceases to be rigorously static and 
only the dynamical picture survives, a picture where the Killing vector with its associated 
energy/surface gravity/temperature plays a much more minor role. 
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hospitality. SAH was supported by the National Natural Science Foundation of China under 
grants 10375081, 10473007 and 10771140, by Shanghai Municipal Education Commission 
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